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Noninvasive, non-contact and m vivo monitoring blood glucose is a long needed pathology tool 
for saving patients from recurring pain and hassle that can accompany conventional blood glucose 
testing methods. Optical coherence tomography known for its high axial resolution imaging modality 
is adopted in this article for monitoring glucose levels in tissue like media non-invasively. Making 
use of changes in reduced scattering coefficient due to the refractive-index mismatch between the 
extracellular fluid and the cellular membranes and armed with a theoretical model, we establish 
a correlation between the glucose concentration and reduced scattering coefficient. The scattering 
coefficients are extracted from the deconvoluted interference signal by using Monte-Carlo simulation 
with valid approximations. A program code using NT LabVIEW"^^^ is developed for automation of 
the experiment, data acquisition and analysis. 



I. INTRODUCTION 

The method of optical coherence tomography (OCT) 
was first introduced by Huang and colleagues. It has 
been widely applied to medical imaging and diagnostics^. 
OCT have the capability to acquire two- and three- 
dimensional tomographic images in biological tissues. 
These applications are limited by penetration depth, 
cross sectional area, dynamic range and signal to noise 
ratio (SNR). However, it is successfully appHed to trans- 
parent ocular organs where light scattering is minimum. 
Multiple scattering, which becomes dominant at large 
depths, is the fundamental limitation preventing OCT 
from reaching a large probing depth in turbid medial'. 
OCT is widely used as a biomedical imaging modality. 
We extent this idea towards a new direction where tomo- 
graphic (cross sectional) imaging is not the prime goal, 
instead we measure the optical properties of stratified 
media with better accuracy. 

The optical properties themselves can potentially pro- 
vide information to monitor tissue metabolic status or 
to diagnose disease. Optical approaches to study tur- 
bid media in the presence of chiral^, components have 
generated interest because of their potential use in non- 
invasive glucose monitoring for diabetes patients. Light 
scattering occurs in tissues because of the mismatch of 
index of refraction between the extra cellular fluid (ECF) 
and the membranes of the cells composing the tissue. In 
the near-infrared region, the index of refraction of the 
ECF is uecf ~ 1.348-1.352 while the index of refraction 
of the cellular membranes and protein aggregates is in 
the range Uceii ~ 1.350-1.460^1^. It is well known that 
adding sugar to water increases the index of refraction of 
the solution. Similarly, adding glucose to blood in turn 
raise the refractive index of the ECF, which will cause 
a change in the scattering characteristic of the tissue as 
a whole. Hence, tissue glucose levels are correlated with 



scattering coefBcients based on changes in the refractive 
index of extra cellular fluid. Of late, measurements on 
light scattering by blood shows promising correlation be- 
tween blood glucose and reduced scattering coefficient^. 
On the other hand, the light scattering is not influenced 
by the red blood cells and other chemical composition of 
blood. At the same time, monitoring of giycemic status 
in patients with diabetes requires determination of blood 
glucose concentration. Significant efforts have been made 
by several groups in the past few decades to develop a 
bio-sensor for noninvasive blood glucose analysis. 

Different optical approaches were proposed to achieve 
this goal. These approaches include polarimetry, Raman 
spectroscopy, near-infrared (NIR) absorption and scat- 
tering, and photo- acoustic o^'^d^d^d^ . Although these 
techniques are promising, they have limitations associ- 
ated with low sensitivity, accuracy and insufficient speci- 
ficity of glucose measurements at physiologically relevant 
levels. Of late, Larin et a l^^d^ have proposed a possible 
change in slope of OCT signal due to changes in glucose 
in blood by optical coherence tomography. However sig- 
nal analysis procedure is quite complex and less accurate. 
They used a used a linear fit model to deduce OCT slope 
by using the least-squares method. We employ a three 
step procedure which is entirely different from their tech- 
nique: (i) OCT signal is deconvoluted from the source 
function, (ii) extraction of optical properties of the tur- 
bid media using Monte Carlo simulations with few valid 
approximations and (iii) interpret the changes in glucose 
concentration from the measured optical parameters. 

Some theoretical models were developed to understand 
the governing physical process and to better interpret the 
OCT signal in highly scattering media. Pan et al, estab- 
lished the relationship between the path-length resolved 
reflectance signal and the OCT signal using linear system 
theor y ^^d^ . Monte Carlo technique was employed to sim- 
ulate the path-length resolved reflectance but could not 
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be able to separate the effects of the singly scattered light 
and the multiply scattered light.— The OCT signal was 
split into (i) summation of singly back-scattered light (co- 
herent) and (ii) multiply scattered light (partially coher- 
ent). The effect of multiple scattering on the formation 
of speckle patterns and the degradation of image con- 
trast were demonstrated. In reality, light scattering in 
turbid media is a complex process, and it is only an ap- 
proximation to assume that the OCT signal is from sin- 
gle back-scattering alone. Photons still contribute to the 
OCT signal after a limited number of scattering events. 
The multiple scattering effects are clearly demonstrated 
in terms of the spreading of the point spread function 
(PSF). 

In the present paper, we employ Optical Coherence To- 
mography to monitor the reduced scattering coefficient 
for different values of glucose concentration in an aque- 
ous solution with Intralipid as the scatterers. Monte- 
Carlo simulation technique with valid approximations is 
adopted to understand the contribution of the multiple- 
scattered light obtainable from OCT. The experimen- 
tal observation are supported by the theoretical analysis 
based upon the transport theory. A strong correlation 
between reduced scattering coefficient and glucose con- 
centration is established. Maximum measurement error 
of 5% is observed at hypo-glyccmic range. 



II. THEORY 



Michclson Interferometer is the basis of any OCT setup 
(as schematically shown in Fig. 1), a motor controlled 
reference mirror and sample with focusing assembly are 
kept in two arms of the interferometer. Light scattered 
from the sample arm and the light reflected from the 
reference arm interfere and is detected by a photodiode. 
The Dopplcr frequency generated by constant scanning 
speed of the mirror is modulated by a coherence func- 
tion of the low coherent light source. A bandpass filter 
centered at the Dopplcr frequency acts as a coherence 
gate for signal detection. The sample in a OCT setup 
used to be a turbid medium such as tissue, plants, com- 
posites, etc. Light scattered from a turbid medium may 
be broadly classified under two categories, viz., (i) the 
least scattered light which undergoes only single or very 
little scattering and (ii) diffusely scattered light, which 
undergoes multiple scattering events. From the princi- 
ple of classical theory of scattering, least scattered light 
maintains coherence whereas multiply scattered light lose 
coherence. Due to finite width of the coherence gat- 
ing employed in OCT, multiple scattered light with path 
length difference falling within the coherence length of 
the source are detected. The optical irradiance at the 
detector is superimposition of all the light fields reflected 
from within the scattering sample and the reference mir- 
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FIG. 1: Schematic of the experimental setup. SLD - Super- 
luminiscent diode, IVI - Mirror, BS - Beam Splitter cube, TS 
- stepper motor controlled translation stage, PD - photodiode, 
S - sample 



ror and is given by. 



/ E',{t,Ls)dLs+Er{t,T) 
Jls 



(1) 



where, r [ = (Ls — Lr)/c = /S.L/c ] is the time delay 
corresponding to the round trip optical path length be- 
tween two beams, Lg[= L^q + X^i^o^s*] takes care of 
round trip path length to the sample surface and the to- 
tal path length within the sample that accumulates dur- 
ing each scattering. E'^{t, Lg) is path length resolved field 
intensity. The first part of cq. ([T]) is important for the 
present study since it contains information regarding op- 
tical properties of turbid media. 

OCT signal obtainable from a turbid media is a convo- 
lution of path-length resolved diffuse reflectance (arises 
from turbidity of the media) and the low-coherence func- 
tion arises from coherence property of laser source. Ac- 
cordingly, OCT signal may be rewritten as. 



h{Lr) = 2/Ur[i?(i.)'/' ® C{L,)]. 



(2) 



Where, Id{Lr) is the OCT signal detected by a photodi- 
ode. Is and Ir are the signals from sample and reference 
arms respectively, R{Ls) is path-length resolved diffuse 
reflectance and C{Lg) is low-coherence function of the 
source. In order to extract the scattered light (R{Ls)) 
from Id, it is deconvoluted with the coherence function 
C{Ls)- The scattering coefficients are obtained by flt- 
ting the deconvoluted signal with the following equation 
obtainable from transport theory 



V^Ud{r) ~ KlUd{r) = -—p^M{r 



47r 



(3) 



Here, Udir) is the average diffuse intensity, Pq is the total 
radiating power and Kd (= yip^J^aPi^r) is a fluid con- 
stant which varies with concentration of the fluid. Here, 
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jjLa, Ms s-iid ^tr are the absorption, scattering and trans- 
port coefRcients, respectively and p is the number density 
of scatterers. 

Light scattering in turbid media depends strongly on 
the value of anisotropy parameter g, which is the average 
cosine of the scattering angle. Total forward scattering 
occurs for <? = 1 while isotropic scattering occurs for g = 
0. In tissue like media the light scattering is generally 
highly forward peaked since the anisotropy parameter 
g > 0.8. Since OCT signal composed of backscattered 
photons, the scattering coefficient is replaced with the 
reduced scattering coefficient 

Ms = (1 - .9)Ms 

instead of the conventional scattering coefficient, to take 
care of changes arising due to the anisotropy in the sys- 
tem. Also, the wavelength of light source can be ap- 
propriately chosen to minimize the light absorption such 
that fi^ « Ms.— Accordingly, the transport coefficient 
is now redefined in terms of the absorption and scattering 
coefficients as 

Mir = Ma + (1 - .9)Ms = Ma + Ms- (4) 

After simplification, the solution of Eq. ([3]) is found to 

be 

--^l^.=^M-Ka{r~r,)]. (5) 
To Udo{ro) 

From the above equation and radial distribution of diffuse 
light relative change of scattering coefficient and glucose 
concentration can be measured. 



III. EXPERIMENT 

We use an IR superluminescent diode (SLD) emitting 
at 810 nm (Hamamatsu) light source. The wave- 
length of the light source is appropriately chosen such 
that the tissue absorption can be ignoredi^. The scat- 
tered light is detected by a photodiode (Hamamatsu) and 
the signal is measured using a lock-in amplifier (EG&G 
7265). 

The motor controlled stages, data acquisition system 
and the lock-in amplifier are controlled through PC. Af- 
ter deconvolution of the OCT signal (J^) with reference 
to the source function {C{Ls)) of fight source, diffuse re- 
flectance signal {R{Ls)) is obtained. It is further used 
with transport theory to calculate the reduced scattering 
coefficient of the turbid media. The data collected for 
different concentration of glucose are analyzed and the 
scattering coefficient is estimated for each measurement 
with valid Monte-Carlo approximations^^. 

Sample description: Glucose solutions prepared 
with different concentration were mixed with turbid me- 
dia (Intralipid 0.1 v/V, values of scattering coefficient 
and anisotropy parameter that are used for the tissue 




FIG. 2: Front-panel of Monte-Carlo simulation program written 
using LabVIEW. The dotted curves are obtained from experimen- 
tal data while the solid lines are fit data using simulation. The 
simulation constants are displayed at the bottom. The results 
are displayed after 1000 iteration (optimized). 

phantom are 50cm^^ and 0.8 respectively), used as sam- 
ple. The exact values of the glucose concentration used 
in experiment ranges from 20mg/dl to 2000mg/dl to es- 
tablish the correlation between the hypoglycemic stage 
(< 64 mg/dl) and hyperglycemia stage (> 100 mg/dl). 
One should note here that the normal blood glucose level 
range is 64mg/dl to lOOmg/dl. The time allotted for 
settlement and interaction of glucose with Intralipid is 
2 minutes. Since, the reduced scattering coefficient of 
glucose are not altered immediately after the addition 
glucose. 

IV. RESULTS AND DISCUSSIONS 

In order to carry out the simulation a program code 
is written using Lab VIE W"^^^ with 1000 iterations (with 
option for N number of iterations). The program window 
is shown in Fig. 2. The program makes use of eqs. ^ - 
([S]) for simulation and data fitting. The method adopted 
by Larin et a\M^ to monitor glucose concentration in 
tissue phantom, obtains tomographic images (2D), which 
are averaged to ID distribution of light in depth. The 
ID distributions were plotted in a logarithmic scale to 
find the slope of the distribution at different depths using 
least-squares method. Which leads a less accurate results 
with large error in measurement. We adopted a different 
direction, where the OCT signal is deconvoluted from 
the source function of the SLD. The deconvoluted signal 
is fitted using MC simulation with the PSF as defined 
earlier. The following parameters are estimated from the 
program code: Msj Ma and fit constants arising due to 
experimental dark noise and background noise. 

The OCT signal obtained for different concentration 
of glucose are exhibited in Fig. [3] for a fixed concentra- 
tion of intralipid. With increase in glucose concentra- 
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tion the amplitude of the OCT signal decreasesJ^ Also, 
one can notice that the Gaussian width of the curve in- 
creases. However, one can not extract much information 
about the scattered light from these curves, since a strong 
coherence function from the low coherent light source 
is convoluted on the scattered signal. The coherence 
data obtained from an OCT signal of the non-scattering 
medium like water is used as the source function {C{Ls)). 




FIG. 3: OCT signal for different glucose concentrations. For 
larger concentrations, the amplitude of signal decreases, while 
the width of the curve increase. The deconvoluted signals are 
shown in inset. The OCT signal obtained for water is taken 
as the source coherence function. The deconvoluted signals are 
normalised. 

Other OCT signals obtained for different concentrations 
of glucose are deconvoluted and shown in the inset of Fig. 
[3l The deconvoluted signals contain more information 
about the light scattering. With increasing concentration 
the area as well as the tail part of the signal increases. 
This is a clear indication of increasing value of reduced 
scattering coefficient with decreasing glucose concentra- 
tion. The measured scattering coefficients obtained after 
Monte-Carlo simulation are depicted as a semilog plot in 
Fig. 31 The curve (solid circles) exhibiting the nature of 
reduced scattering coefficient for different concentrations 
of glucose demonstrates a near logarithmic nature. In 
hypo-glycemic region it shows a sharp change while in the 
hyper-glycemic region it has smaller slope leading to less 
accuracy. A similar nature is also observed while mea- 
suring the area of the deconvoluted signal (open circles). 
This also experiences a near logarithmic behavior. Re- 
peated measurement shows that the slope of the curves is 
a constant. Also, addition of more or less intralipid does 



not change this behavior. However, in order to find the 
value of unknown concentration of glucose, one needs an 
initial value at particular concentration of glucose. Since, 
the slope remains a constant, this method could be better 
solution for non-invasive, non-contact, in- vivo monitoring 
of blood glucose concentration. In order to realize this 
technique as a clinical tool more efforts are required. 
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FIG. 4: Semilog plot of reduced scattering coefficient and curve 
area with glucose concentration. These data extracted from 
Monte-Carlo simulation and fitting it to the experimental data. 
The error bars are obtained after 25 measurements. 

V. CONCLUSION 

To conclude, armed with light scattering technique and 
optical coherence tomography, we make an attempt to 
study a correlation between, glucose concentration and 
reduced light scattering coefficient which is better and 
accurate with than Larin et a^^^. A semilog plot of 
glucose concentration with reduced scattering coefficient 
suggests a linear relation. The value of reduced scatter- 
ing coefficients are determined using Monte-Carlo Sim- 
ulation. The technique promises to be a possible clini- 
cal tool for non-invasive and non-contact measurement of 
blood glucose. However, to predict exact value of blood 
glucose without feeding any initial parameter, the equip- 
ment needs large number of training in a pathology lab- 
oratory. 
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